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Abstract. Finite thermal velocity modifications to electron slowing-down rates may be important for the deduction of solar flare 
total electron energy. Here we treat both slowing-down and velocity diffusion of electrons in the corona at flare temperatures, 
for the case of a simple, spatially homogeneous source. Including velocity diffusion yields a consistent treatment of both 
'accelerated' and 'thermal' electrons. It also emphasises that one may not invoke finite thermal velocity target effects on electron 
lifetimes without simultaneously treating the contribution to the observed X-ray spectrum from thermal electrons. We present 
model calculations of the X-ray spectra resulting from injection of a power-law energy distribution of electrons into a source 
with finite temperature. Reducing the power-law distribution low-energy cutoff to lower and lower energies only increases the 
relative magnitude of the thermal component of the spectrum, because the lowest energy electrons simply join the background 
thermal distribution. Acceptable fits to RHESSI flare data are obtained using this model. These also demonstrate, however, that 
observed spectra may in consequence be acceptably consistent with rather a wide range of injected electron parameters. 
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• *"H 1. Introduction 

;_i ■ X- and y-ray radiations give the most direct window on accelerated electrons in flares. They have revealed that accelerated 
5^ ! partic les, electrons and/or ions, are an energetically major product of the flare energy release process (e.g. lVilmer & MacKinnon! 
I2003I). 

iBrown et alJ J2003I) have emphasised the importance of the source-averaged electron distribution as a useful 'halfway house' 
between the observed photon spectrum and the distribution of electrons initially injected into the source region, i.e. the immediate 
product of the acceleration process. Assumptions about the dominant factors in electron transport then allow deduction from 
the source-averaged electron distribution of the distribution output by the acceleration process. Quantities like the total energy 
rel eased in the f orm of fast electrons follow immediately. 

iBrownl dl97ll) first analysed the case in which electrons slow down via Coulomb collisions in a cold target, i.e. a region in 
which ambient particle thermal speeds are all very much less than those of the X-ray emitting electrons. Key results were given 
for a photon spectrum 7(e) (photons cm~ 2 keV -1 s _1 ) depending on photon energy e as a power-law, i.e. 1(e) ~ e~ y for some y > 0. 
Such a photon spectrum is often observed. Assuming a cold target from which no electrons escape, it implies an injected electron 
energy distribution depending on electron energy E as E~ r ~ l . The total energy content of such a distribution is governed by the 
lowest electron energy for which this power-law form holds good. Unfortunately observations remain ambiguous on the likely 
value of this lower cutoff, so the total flare energy in accelerated electrons remains uncertain by more than a n order of m agnitude. 
The flare energy in electrons of energies > 25 k eV appears to be a large fraction of the total flare energy dLin & Hudsonlll976l 
iHovng et al.lll976l Isaint-Hilaire & Benzl Eo05): observations even exist suggesting a low energy cutoff in the 2-5 keV range 
jKaneet all 19921) 
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Emslie] J2003I) has pointed ou t that the cold targe t assumption may be invalid for the lowest energy (few keV) accelerated 
flare electrons. Spatial structure (Emsl ie et al.l l2003) of RHESSI (Reuven Ramaty High Energy Solar Spectroscopic Imager) 
images in particular suggests that these electrons stop entirely in the corona, in high temperature (> 10 7 K) regions. Then the 
test particle slowing-down rate no longer increases monotonically with decreasing particle energy; as particle speed approaches 
ambient particle thermal speeds from above, the rate of loss of energy to background particles decreases, exhibiting a zero for 
a fast particle energy E cr n very close to the ambient electron thermal energy. Emslie contends that this value should be used as 
the minimum possible lower energy cutoff when evaluating fast electron total injected energy; electrons below this energy do 
not slow do wn monotonically , instead merging with the background thermal distribution. Emslie's suggested procedure has been 
employed bv lLin et"al] ll2003h to estimate the fast electron energy content in the flare of 23rd July 2002. 

Emslie's discussion is expressed entirely in terms of the systematic slowing-down rate. This gives valuable insight but cannot 
give a complete description of the evolution of injected electrons. Suppose we inject a mono-energetic electron distribution at 
Zscrit- Clearly, although the systematic slowing-down rate at E cr \t is zero, the injected electrons will not stay indefinitely at E cr \ t ; 
they will spread out in energy in such a way as to eventually join the ambient Maxwellian population, doing so in the first instance 
via diffusion in velocity rather than systematic slowing down. A more complete treatm ent is needed to discuss the form of the 
photon spectrum and what it is telling us about flare fast electrons. iMcClementsI i 19871) has included velocity diffusion effects, 
but only as one component of a complicated treatment which also features a number of other processes, and he does not explore 
the issues we address. 

In this contribution we examine the evolution of injected electrons when the cold target assumption breaks down, in the 
slightly idealised case of a homogeneous source and constant background temperature. We include velocity diffusion as well as 
systematic slowing down, and reformulate the interpretation of observed photon spectra. The next Section formulates the problem 
and gives some analytical discussion. Section [3] illustrates our discussion with some calculated spectra, compared to RHESSI 
data. Section@]gives brief conclusions. 

2. Velocity diffusion of fast electrons 

2.1. Assumptions; Fokker-Planck equation 

In order to illustrate the consequences of velocity diffusion for photon spectra we consider an idealised problem in which all 
injected electrons thermalise in a uniform, homogeneous medium, characterised by a single, ambient electron density n e and 
temperature T e . Loss at boundaries will have a negligible influence on the electron distribution function and pitch-angle in- 
formation will not be important for the calculation of the total, emergent photon spectrum, particularly since bremsstrahlung 
directionality is unimportant at the few keV photon energies appropriate here. We can gain significant insight, and also solve 
a problem appropriate to understanding the X-ray emission integrated over the whole of the event, by studying a steady-state 
situation, so no quantity depends on time. In practice n e and T e will evolve as a result of the thermal and hydrodynamic response 
of the atmosphere to the flare energy release, but this is a complication of detail rather than principle and we ignore it in the 
interests of gaining insight. (Relevant electron timescales such as the thermalisation time are at most on the order of seconds, 
whereas bulk changes to the plasma, such as changes of temperature as characterised by variation of the soft X-ray flux, take 
place over timescales of a few minutes.) Thus we can characterise the (pitch-angle integrated) electron distribution everywhere 
in the source by a single function /(v) ((cm s -1 ) -3 ) of velocity v (cm s _1 ). The normalisation of / is given by 

/-»oo 

An f(v)v 2 dv = N e , (1) 
Jo 

where N e is the total number of all electrons in the source. 

An electron of 10 keV initial energy will stop in a column depth of 2 x 10 19 cm 2 of fully ionised hydrogen, inside the 
coronal portion of a loop, e.g. for densities > 2 x 10 10 cm 3 and loop lengths > 10 9 cm, conditions not infrequently inferred in 
flares. Thermalisation, in the alternative case that this energy is close to the thermal speed, will occur in a comparable distance. 
In addition, magnetic field convergence may further enhance the coronal residence time of electrons and increase the effective 
distance available for thermalisation. Electrons well above thermal speeds will experience cold target conditions throughout the 
corona and chromosphere and will in any case be described correctly by what follows. Thus, while possibly not the case in all 
events, it is not unreasonable that all electrons for which finite thermal velocity ('warm') target effects are important thermalise 
in the coronal, warm target region. 

This steady-state treatment will be valid as long as we use it on times that are not long enough for the injected electrons 
to become significant in number compared to the thermal distribution, but are longer than relaxation times for most injected 
electrons. Alternatively, we may regard it as giving the time integral of the distribution functi on in the case of a n initial, impulsive 
injection, in which case the source function S is actually the initial condition on / (MacKinnon & Craig 1991). The time integral 
of / is the necessary quantity for calculation of the total bremsstrahlung photon yield. 

We use the Fokker-Plan ck formalism for particle evolution under binary collisions (e.g. iRosenbluth et alJ Il957l 
llvlontgomerv & Tidmanl 19641) . We also make the assumption that the fast particles are 'dilute', in the particular sense that they 
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may be ignored in calculating the velocity space drift and diffusion coefficients: these m ay be ev aluated purely from the back- 
ground distribution. Then the steady-state, Fokker-Planck equation for f(v), derived from Helander & Sigmar (2002, p.37-38), 
may be written as 

' C> ^(v)-v$'(v)}(/4?)) =S(y), (2) 



v 2 dv \ ( 2v dv 

where the function <t> is the error function: 

2 C x 2 

<D(x) = erf(jc) = — e~ r dy. (3) 
yn Jo 

Here velocities v have been normalised to the ambient electron thermal speed vj - y/2kT/m e , and times to the electron 
thermal collision time t c = 47re/jm 2 v^«~ 1 e~ 4 (ln A) -1 Although there is no time-dependence in this problem, the electron injection 
function S (v) is of course per unit time. 

Equation (0 needs two boundary conditions. The boundary condition at infinity is: (/ + l/2v df /dv) — > as v — > oo. This 
condition ensures that there is no flux of particles out of the system at infinity. For the other boundary condition we fix f(v) at 
v = 0: /(0) = /o, consistent with the conditions for validity of our steady-state treatment and with the assumption of 'diluteness' 
that justified the linearisation of the Fokker-Planck equation, /o describes the background thermal distribution. We integrate 
Eq. once from v to infinity and use the boundary condition at infinity, then again from to v, employing an integrating factor 
e v and using the boundary condition at 0. Thus we find the solution 

; (>' 2 -' 2 ) 

(4) 



/(v) =/ e- +2 - d/ v" 2 S(v")dv", 

Jo O(v') - v'O'O') J,, 



which will be used in Sect.|3]to calculate distributions and resultant photon spectra for various forms of S . The result|4]is rendered 
rather impenetrable, however, by the function <J>. An illuminating, semi-quantitative analytical discussion may be carried out by 
invoking the large argument form of O(v), strictly applicable for v » 1, in which <£(v) approaches unity for large v. In this case 
the Fokker-Planck equation becomes: 

v 2 dv if 2v dv) ^ ^ 
2.2. Approximate solution 

First, we note that the LHS of Eq. may be rewritten 



2v 2 J dv 2v 3 dv 2 ' 

showing explicitly that the systematic rate of change of v (the coefficient of df/dv) does i ndeed display a zero, even in this 
approximate form, at v = 1/ V2, quite close to the zero found using the full form of <t> by lEmsliel d2003l) . All the necessary 
qualitative features are included in the description of Eq. (|5}, in spite of its approximate nature, and its solutions will have the 
appropriate qualitative properties. 

Note that the warm target corrections to the sy stematic slowing-down rate and the dispersive term both become important 
in the limit v — * 1/ y2. The arguments of Emslie (2003) rest on the presence of the zero (here at v — 1/ y2) in the systematic 
slowing-down term. Electron slowing-down times approach oo, so a finite emergent photon spectrum demands S — > as v — > 
1/ V2. However, in this limit the dispersive term has become important, removing the divergence in electron 'lifetime'. Using the 
boundary conditions as before, and changing the order of integration in the resulting integral, Eq. © has the solution 



Jov poo 
u 2 S(u) (e" 2 - l) du + (l - e~" 2 ) I u 2 S(u 
Jv 



)du. (6) 

In the absence of any source S, Eq. (0 has the background Maxwell-Boltzmann distribution (/oe~ 1,_ ) as its solution, as 
indeed does Eq. (|2)- The description in terms only of systematic slowing-down rates divorces the fast particle and background 
distributions. This is no longer the case in this diffusive treatment: the presence of the boundary condition at v = 0, which must 
be satisfied using the same background density used to calculate the drift and diffusion coefficients, ensures that the fast particle 
distribution merges smoothly with the thermal 'core'. It follows that we are obliged to include also the contribution to photon 
emission from the thermal plasma, if we are indeed looking at photon energies such that velocity diffusion is important for the 
emitting electrons. 
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In the limit v — > oo, Eq. (|6j becomes 

u z S(u)du. (7) 

Recall that our f is identical with the mean source electron distribution, the key quantity in interpreting X-ray emission 
jBrown et al.l EoQ3). With the traditional assumptions of fast electrons slowing down in a cold, thick target, this mean distri- 
bution is just the cumulative distribution of the injected energy distribution. Eq. |7} reproduces this result, as it should. We can 
rapidly recover well-known results in that limit, for instance Brown's (1971) relations between energy power-law spectral indices 
of observed photons and injected electrons. 

We see that the three terms in the solution consist of: the Maxwell-Boltzmann core of the distribution; a term which 
resembles the cold target result more and more closely as v — > oo; and a term which forces these two regimes to merge smoothly. 

2.3. Mono-energetic injected population 

The special case of a mono-energetic form for S is instructive: 
S 6(v - Vg) 

S(y) = , (8) 

for some velocity v . Then the solution becomes: 

ff/o+Sofe-'o-l^e-" 2 ifv>v , 

/(v) = I V c (, U f (9) 
[foe +S [l-e 1 J if v < v . 

For v < vo the distribution is composed of the original, background Maxwellian distribution plus a component which is iden- 
tical to the cold target result for v, vo » 1, but which approaches as v — > 0. This additional, non-Maxwellian component 
becomes less and less significant for smaller and smaller vo- For v > Vo, the distribution is identical with the original background 
Maxwell-Boltzmann distribution, only with its normalisation increased. Since we must have So <z fi> for validity of the original 
linearisation of the Fokker-Planck equation, we see that the distribution will resemble the original Maxwellian more and more 
closely as vo gets closer and closer to 0. This justifies the qualitative comments made in Sect. ^ electrons injected close to the 
thermal speed diffuse in energy rather than slowing down monotonically, merely adding their number to the original background 
Maxwellian. Figure^iHustrates this. 

2.4. Deduction ofS for a power-law photon spectrum 

As mentioned in Sect. HI the photon spectrum may in principle be inverted to yield the mean (source-averaged) electron distribu- 
tion jBrownll 197 It Brown et alJl2003l) . This is identical with the distribution function / in the special case of our homogeneous 
source. If observations have given us / in this way, Eqs. (0 or Q then immediately give S . Consider the case of a power-law 
photon spectr um 1(e) ~ e~~> ' . Assume for the moment that this form holds at all photon energies of interest. Then the results of 
Brow5 lll97ll) give us /(v) ~ v~ 2r ~ 3 . Inserting this form into Eq. (|5jl we find 

(10) 



S(v)~(2 r+ 3)v-^{l-^i 



As in the case retaining only systematic energy loss (Ems liel2003l) . S has a zero, changing sign at 



Vr + 5/2. 



S(v) takes negative values for v < v*. We have seen that a diffusive treatment underlines the necessity of including the radiation 
from the thermal, 'core' part of the distribution. In assuming the power-law photon spectrum to be appropriate at all photon 
energies we have implicitly neglected this contribution, v* represents the lowest velocity at which the single, uninterrupted 
power-law photon spectrum can still be reconciled with the presence of the Maxwell-Boltzmann core. Below v„ we would have 
to actually remove particles from this core to prevent deviations from a power-law photon spectrum; hence v„'s dependence on 
y. Moreover, we cannot 'overcome' the core Maxwellian distribution by, for example, injecting a power-law energy distribution 
of electrons that persists down in energy towards thermal speeds. As we saw in Sect. 12.31 these electrons mostly thermalise 
diffusively, producin g only a slight modification to the core. 

We might follow lEmsliel d2003l) and evaluate total electron energy content by integrating S , given by Eq. dlO> . from v* to 
oo. Rather than adopting a lower energy cutoff for the power-law which evidently holds at high energies, however, this approach 
underlines the need for a consistent treatment of radiation from both thermal and accelerated electrons. 
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Fig. 1. Form of the relaxed, steady-state distribution resulting from injected, monoenergetic electrons. The different curves repre- 
sent different values of the velocity (in units of the thermal speed) of the injected electrons. The noticeable change in behaviour 
between the curves for vo = 2 and vo = 3 highlights the injected electrons' transition from being in the thermalisation regime to 
being properly non-thermal. 

3. Numerical illustrations and application to data 

3.1. Numerical illustrations 

We return now to the full solution of the Fokker-Planck equation as given in Eq. ©, and provide some illustrative examples 
relevant to solar observations. We adopt as the source function S (v) a power-law, employing the form 

S(v) = S(v)H(v-v ) 

= S (6 v -3)v s ( ;- 3 v- s 'H(v-v Q ), (11) 

where § is normalised such that per unit time (normalised to the electron thermal collision time) there are Sq particles injected 
in total at velocities above Vo, and S is prevented from going to infinity at low velocities by Heaviside's step function H which 
removes all particles with velocities less than vo. 

For a homogeneous source, the emission rate of photons of energy e per unit energy range per unit volume, dj/ de (photons 
s" 1 keV" 1 cm" 3 ), may be found by multiplying the distribution function by vdcr/de to obtain the instantaneous rate of emis- 
sion of photons by electrons in the velocity range v — > v + d 3 v, then integrating over all velocities ( Brown| ll97ll) . noting that 
d 3 v = 47rv 2 dv. This gives 

dj , f°° do- , 

-f =n e v\ /(v) — v 3 dv, (12) 

where n e is the background plasma number density and ^ is the Bethe-Heitler cross-section: 

dcr _ Q m e c 2 1 1+ VI -e/E \ 

de eE \\ - yJT^ljE) ' 

Here, E is the electron kinetic energy and 2o is given by 

8 2 
Qo = jar e , 
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Fig. 2. Simulated X-ray photon spectra, showing the effect of altering the ratio of fo to Sq. The photon energy is normalised to 
the kinetic energy of an electron of thermal velocity. The emitted intensity values are arbitrary, since we are interested only in 
relative changes to the spectral profiles. 



where the fine structure constant a w 1/137 and r e = 2.82 x 10 13 cm is the classical electron radius. The photon spectrum, 
d// de (photons s _1 crrT 2 keV -1 ), that would be observed by RHESSI is given by 



de 



V dj 
Anri de 



(13) 



where V is the volume of the source and r m is the distance from the Sun to the Earth. Since radiation from the whole emitting 
volume is observed, the value of V will be determined implicitly by the spectral fitting process (see Sect. 13 .2I > and need not be 
separately evaluated. 

As previously stated, the prese nce of O i n Eq. I@J renders a full analytical solution intractable. Therefore we proceed numer- 
ically, using Romberg integration ( Pres s et al.ll992l) to evaluate Eq. J12I with / given by Eq. (0J. 

The parameters which may be varied in the numerical simulations are: the ratio of fo to So i-e. the relative magnitudes of the 
background Maxwellian population and the injected power-law electrons; the lower cutoff velocity of the injected electrons, vo; 
and the spectral index 6 of the power-law. It should be noted that in the majority of the literature, the power-law used to model 
flare electrons is a power-law in energy of the form S (E) = S qE~ 6 . For consistency and ease of comparison we shall refer to this 6 
in subsequent discussion. The corresponding 8 V for our velocity power-law, Eq. (II 1> . is related to 6 by the expression 5 V -26-1. 
Unless otherwise stated, the default values of the parameters are fo/So = 10 8 , 6 = 4.0, and vo = vr. 

Figure 13 illustrates the effect of altering the ratio fo/So- The logarithmically-plotted photon spectra consist of two main 
regions: a straight power-law profile at high photon energy, blending smoothly into a Maxwellian profile at lower energy. As 
would be expected, increasing the relative contribution of the Maxwellian background has no effect at high photon energies since 
here the profile only contains contributions from electrons of the photon energy or higher. However, a larger fo/So value leads to 
a correspondingly higher contribution to the Maxwellian portion of the spectra from the background plasma. Furthermore, this 
larger value also corresponds to an increase in the photon energy up to which the Maxwellian impinges on the otherwise straight 
power-law: the profile departs from the straight portion at higher energy for a larger fo/So- 
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Fig. 3. As Fig. [2 but with variation of the power-law spectral index 6, and over an expanded range of photon energy. 



The alteration of the electron energy spectral index 5 is depicted in Fig. [3] As may be seen, this has minimal effect at low 
photon energy, but a larger 5 results in a correspondingly steeper slope in the power-law region of the spectrum. A greater 
value for 5 also causes a more rapid reduction in the total number of electrons as a function of increasing energy in the injected 
population. Thus, a higher 6 leads to a relative reduction in the intensity of the power-law spectrum for a given photon energy. 
This reduced contribution from the power-law electrons also increases the photon energy up to which the Maxwellian element 
forms a significant part of the resultant profile. Consequently, the departure from the straight power-law profile occurs at a higher 
photon energy for larger 5 values, visually mimicking a non-existent change in the power-law low energy cutoff. 

Actual variation of the cutoff, vo, is shown in Fig. |4] Since electrons injected below a few times Vj thermalise rather than 
slowing down systematically, allowing the cutoff to extend to lower energies merely adds electrons to the 'background' Maxwell- 
Boltzmann distribution. This explains the counterintuitive result, clearly visible in Fig. @J that a lower value of vo results in 
a spectrum which attains power-law form at higher photon energies: the large number of injected electrons at low energies 
thermalise and enhance the Maxwell-Boltzmann distribution, concealing the lower-energy portion of the power-law form. 



3.2. Comparisons to RHESSI data 

The recent launch of the Reuven Ramaty High E nergy So lar Spectroscopic Imager (RHESSI) has opened a new era in high reso- 
lution X-ray spectroscopy of solar flares dLin et all2002h . T he analysis of solar flare spectra has revealed statistically significant 
deviations from a simple isothermal and power-law mo del dK ontar et afll2003l) . The detailed analysis of the X-ray producing 
spectra using a model-independent inversion technique dPiana et alJ 20031) shows deviations from the pure isothermal model. 
These new findings can be treated as the manifestation of velocity space diffusion in warm target plasma. 

For illustrative purposes, we consider a few example events with sufficiently high count rate to provide reliable photon 
statistics. We have limited our analysis to within the energy range 10-50 keV, where thermal and nonthermal components merge. 
Below 9 keV, the bremsstrahlung continuum is contaminate d by a comple x of strong iron lines. Above ~ 50 keV, spectral features 
not related to the model discussed become dominant ( Kon taret al.l2 003). 
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Fig. 4. As Fig. [2 but with variation of the power-law cutoff vq. 



Table 1. Optimal fit parameters for our model fits to RHESSI data. 



Flare 


A/So 


vo 


T e (MK) 


6 


Efast (J) 


£therm (J) 


E lot (J) 


2002/03/17 


1.46xl0 7 


1.18 


21.2 


6.30 


3.79xl0 25 


1.43xl0 24 


3.94xl0 25 


2002/08/21 


3.40xl0 7 


2.66 


29.7 


4.86 


1.47xl0 23 


1.38X10 23 


2.84xl0 23 


2002/08/22 


5.28xl0 6 


2.90 


25.0 


5.69 


3.75xl0 23 


5.32xl0 23 


9.07xl0 23 



fifast and Siherm are the energy contents of the fast and background electrons respectively, calculated using the optimal parameters given. 
E iot is the total energy content of all the electrons. 



We fit the model spectra to the observed spectra by optimising over 4 parameters: the relative magnitudes of the background 
and injected electron populations, fa /So', the injected electron power-law low-energy cutoff, v ; the background temperature, T e ; 
and the injected electron power-law spectral index, 6. The optimisation seeks to minimise an un-normalised^ 2 fit statistic - in 
this case absolute values of x 1 must be treated with caution since the process of deconvolving the RHESSI instrument response 
from the observed counts spectrum to produc e the photon spectr um introduces an element of error on each photon spectrum data 
point which is difficult to quantify precisely JSmithet aljE 002). However, we only compare relative values of the fit statistic to 
optimise the model fits, so un-normalised x 2 is sufficient for our purpose. As may be seen from Fig.|5j sets of optimal model 
parameters may be found which give model spectra that closely match the observed RHESSI spectra. The model parameters 
corresponding to the smallest fit statistic in each case are given in Tabled 

As previously discussed, our steady-state treatment is valid on timescales longer than the collisional timescale but shorter 
than the timescale for changes in the temperature of the flaring loop. Thus, the fitted model parameters describe the flare plasma 
at any instant, but will change with time as the plasma evolves during the flare. To obtain a simple estimate of the energy content 
of the electrons in the example flares, we fit a model spectrum to an observed RHESSI spectrum from during the impulsive 
phase, and multiply the instantaneous energy content by the duration of the phase. To obtain the instantaneous energy, we insert 
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Fig. 5. X-ray photon spectra recorded by RHESSI during the flares of 17th March 2002 at 19:30:48 UT, 21st August 2002 at 
01:38:28 UT, and 22nd August 2002 at 01:50:32 UT. Model fits to the observed spectra are shown as solid lines. 
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Table 2. Optimal fit parameters for fits to RHESSI data using a simple thermal plus power-law model. 



Flare f /S 



T e (MK) 



S Efast (J) Efherm (J) £" t ot (J) 



2002/03/17 2.52x10 s 
2002/08/21 3.98x10 s 
2002/08/22 l.llxlO 7 



3.10 
1.93 
3.16 



22.0 
29.8 
25.2 



6.41 2.68xl0 21 1.44xl0 24 1.45xl0 : 
4.98 8.62xl0 21 1.53xl0 23 1.61xl0 : 
5.91 6.73xl0 22 5.37xl0 23 6.05xl0 : 



£f asl and -Btherm are the energy contents of the fast and background electrons respectively, calculated using the optimal parameters given. 
£tot is the total energy content of all the electrons. 



the optimal fitted values of the relevant parameters into the source function, Eq. dl It . and calculate the total energy content of 
the fast electrons by integrating the electron kinetic energy over all possible electron velocities: 



The electron energy content of the thermal background plasma may also be obtained from the optimal fitted parameters: 



(For our present illustrative purposes, we assume a background plasma density of 10 15 m -3 - typical of the lower corona - to 
obtain the absolute value of /o from the value of the emission measure as determined from the fits.) The total energy in all the 
electrons is the sum of Efo St and E t h erm . 

We find single values of the source parameters to represent the whole of the data time interval. The main assumption in doing 
this is that background parameters (Jo, T e ) do not change. If this is the case then interpreting the time integral of the data gives 
us the same result as integrating a temporal sequence of data fits (as mentioned in Section [2J see also MacKinnon and Craig 
1991). Although these parameters may change, this will be partly because of the relaxation of the fast electrons. Qualitatively, 
the procedure here may overestimate the injected electron distribution, because temperature will increase, and thus more of the 
observed photon spectrum will be due to 'thermal' electrons, as time goes on. To address this issue quantitatively we would have 
to drop the linearisation of the Fokker-Planck equation, Eq. 0, resulting in a considerably more complex problem that we do 
not address here. Such a fuller treatment would also allow us to precisely determine the realm of validity of our linearisation. 

The time intervals we use include the bulk of the hard X-ray emission from the flares in question, but of course a more 
complete discussion of flare energetics would also integrate over the entire history of the flare. 

Also given in Tableware the total energy contents of all electrons and of the injected and background electron populations 
individually for each event studied. The flares of 21st and 22nd August 2002 have comparable energies in the thermal and fast 
components, and a total energy consistent with an M class flare. Both also have low-energy cutoffs at only a few times the thermal 
speed, indicating that velocity diffusion will be relevant. The optimal fitted parameters for the flare of 17th March 2002 result 
in a much larger total energy which would correspond to a larger and more energetic flare. However, in this case the calculated 
energy content of the injected electrons is many times greater than that of the background thermal electrons. This arises because 
the fitted low-energy cutoff is very close to the thermal speed itself, and the power-law spectral index is very large, meaning that 
the injected population will have a huge number of electrons very near to the thermal speed. These will rapidly thermalise and 
give rise to the bulk of the Maxwellian portion of the photon spectrum as shown in Fig. [5] dominating the emission from the 
background plasma. While this set of model parameters corresponds to a minimum fit statistic, and produces a model spectrum 
which closely reproduces the observations, they also imply a situation where the injected electrons are no longer 'dilute' and our 
linearisat ion i s no lo nger applicable. 

iLin et alJ J2003) employed Emslie's (2003) formulae for flare electron energy content in their analysis of the 23rd July 2002 
flare. This analysis assumed that the injected electrons had a power-law low-energy cutoff at approximately the thermal speed 
(T x 23 MK) and, unlike our treatment, does not account for any thermalisation of these lower energy injected electrons. Lin et 
al. found the fast electron energy content for this X4.8 flare to be in excess of 10 27 J. This is very much greater than the highest 
total energies ever deduced for the largest flares. We also predict very large energies in the injected electrons for cases where Vo is 
very low. However, our estimates are not as extremely high as those made using Emslie's formulae, since our treatment includes 
the appropriate velocity diffusion effects for the lower energy injected electrons. 

Due to the unambiguous nature of the straight portion of the spectral profile at high energies, the fitted value of 6 is well- 
constrained. However, as is evident from Figs.[5]andE] variation of the values of fa/So and Vo lead to similar variations in the 
shape of the resulting spectral profile. This suggests the possibility of a degeneracy in the fitted values of these parameters, which 
is indeed the case. The quoted value of the fitted Vo given in the table is that for the most optimal fit. However, it was found that 
the value of Vo could be varied from around one half to two times the optimal value for only a 10% increase in the value of the 
fit statistic. Therefore we are reluctantly forced to conclude that the value of vo is less well-constrained by the data than Emslie's 




(14) 




(15) 
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original argument might suggest. Similarly, fo/So, and to a lesser extent T e , cannot be unambiguously determined. The value for 
the total electron energy varies by around an order of magnitude when the value of Vo is varied over our selected 10% range of fit 
statistic acceptability, with the total energy decreasing as vo is increased. Thus, while the model can reproduce observed photon 
spectra, its nature may preclude a precise determination of the flare electron energy content. 

For comparative purposes, Table [2 gives fitting parameters as derived from fits using a 'simple' t hermal plus power-law 
model, as may be employed in e.g. the OSPEX package in the standard RHESSI analysis software (e.g. lSchwartz et al1l2002h . 
As measured by our fit statistic, these fits are statistically acceptable at a similar confidence level to our model fits. As may be 
seen, the simple fits give consistently larger values of fo/So (and comparable or slightly greater temperatures) than our model 
fits, since all of the thermal part of the spectrum must be accounted for by the Maxwellian background with no contribution 
from thermalising fast electrons. This results in higher energy contents for the thermal background electrons, but the total energy 
contents are lower than for our fits since the thermalisation process is an energetically expensive way of producing 'thermal' 
photons. The simple fit values for the low-energy cutoffs are comparable with those of our model in that they are a few times the 
thermal speed, but again this parameter is not well-constrained. The largest difference is for the 17th March flare, for which our 
model suggested an extremely low vo. The simple fit is more conservative, and consequently does not give the very large energy 
content in the injected electrons as derived from our model fits. 

4. Conclusions 

We have seen that for the case of accelerated flare electrons impinging on a warm target, the effects of velocity diffusion should 
not be neglected. The process by which the lower energy electrons of the injected population thermalise and merge with the 
ambient thermal background is of importance when describing the behaviour of the electrons in this regime. In particular, it has 
the effect of 'smearing out' this region of the resulting bremsstrahlung photon spectrum, which in many cases is therefore not 
well described by a simple isothermal and power-law model. However, this region can be modelled effectively by a treatment 
including velocity diffusion effects. 

A consequence of including velocity diffusion in the analysis is that simple interpretation of observed photon spectra can 
be deceptive. For example, as has been shown in Figs.[2E]and|3] determining the photon energy down to which the spectrum 
remains power-law-like does not allow a simple evaluation of the parameters of the injected electron population, most particularly 
its low energy cutoff. This in turn hinders determination of the flare electron energy content. 

The surprising behaviour of the photon spectra is highlighted by the fits to RHESSI data in Fig. [5] visual inspection of the 
spectra would seem to suggest that the 17th March flare should have the highest temperature background plasma, since the 
Maxwellian portion of the spectrum is more prominent and extends to higher photon energy in this flare than in the others con- 
sidered. However, the fitted parameters imply that the 17th March flare actually has the lowest background plasma temperature. 
The form and extent of the Maxwellian component of the photon spectrum in this case is completely dominated by thermalising 
electrons from the lower energies of the injected population. This result emphasises that the thermalising process effectively 
couples the injected and the background populations, such that their contributions to the overall photon spectrum cannot easily 
be separated. In effect, the distinction between the background and injected populations becomes rather arbitrary at these low 
energies, and it is no longer meaningful to distinguish between a background thermal electron and an electron which has ther- 
malised out of the injected population. Further, because of this strong coupling it is also not possible to 'swamp' the thermal 
region of the photon spectrum by contriving a large injected power-law population with a very low cutoff energy, as the lowest 
energy electrons will inevitably thermalise. However, this thermalisation process is an energetically expensive way to produce 
thermal emission. 

The steady-state solution presented here has illuminated many of the consequences of velocity diffusion in the context of 
solar flares. However, an approach which explicitly accounts for time-dependence would be an interesting further development. 
This would allow modelling of the evolution of the plasma parameters over the duration of a flare, and may also be of benefit for 
more precise evaluations of the flare energy budget. 
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